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Abstract 


As analog and mixed-signal VLSI designs are becoming popular, an accurate mod- 
eling of the subthreshold region for short channel MOSFETs is becoming extremely im- 
portant. This region of operation is exploited in modern application areas, such as space 
electronics, laptops, medical applications, communications, multimedia, etc. Operation of 
the devices in their subthreshold region reduces the total power dissipation, which is crucial 
as currently the packing density is increased manyfold due to the rapid advancements in 
the photolithography techniques. Experimental data on short channel MOSFETs reported 
in literature show that the subthreshold slope is a function of the applied drain voltage, 
which previous authors have failed to model. In this work, we have taken up an existing 
subthreshold slope model (which is independent of the drain voltage) and have modified it 
to include the effect of the drain voltage. Additionally, our model accounts for the effects of 
the effective channel length and the body voltage on the subthreshold slope. We have also 
attempted to represent the fudge factor, used widely in the expression for the characteristic 
length, by a more physical representation than done before. The subthreshold slope model 
has been put in an existing drain current model, valid for all regions of operation, in order 
to obtain the DC characteristics. The simulated results are compared with those reported 
experimentally for MOSFETs having submicron channel lengths, with the minimum being 
as small as 0.075 /xm. The results show an excellent match between the two. 
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Chapter 1 


Introduction 


Invention of the transistor launched a new era of electronic sciences, which created an 
enormous impact on the modern scientific applications. Very Large Scale Integration 
(VLSI) is the assembly of millions of such transistors on a single chip known as the 
Integrated Circuit (IC). These transistors work flawlessly to realize the modern elec- 
tronics applications in multifarious areas such as multimedia, communication, defense, 
etc. The advancements in the VLSI technology over the past two decades have sig- 
nificantly influenced every aspect of electronic design, computers and communication 
systems. 

One of the reasons for the popularity of the VLSI systems is the advancements 
in the fabrication technology which resulted in a very high packing density by reducing 
the size of the individual transistors on the VLSI chip. The MOS VLSI has always 
been dominating in the VLSI application area because of continuous reduction in the 
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size of the MOS transistor. The immediate outcome of the very high packing density 
is the tremendous increase in the power dissipation in the VLSI chips which needs to 
be kept within limit. The emergence of products (ICs) having device count in excess of 
5 millions using CMOS and BiCMOS technologies and having processing power of the 
order of 100 millions operations per second or more are challenged by the high-power 
dissipation. This ultimately led the scientists towards the low power design alternative. 
Now in the current VLSI scenario, the high speed/low power chips are most viable. 
The examples of low voltage VLSI systems are wireless communication chips, memory 
chips and many more battery operated equipment. 

The MOSFET (Metal Oxide Semiconductor Field Effect Transistor) is no 
longer a device suited only for the digital applications. In the past few years, the 
semiconductor technology has witnessed the appearance of many MOS chips combin- 
ing both digital and analog functions. Such ICs are known as the mixed-signal ICs. 
This trend is likely to continue, since most signal in need of processing are analog at the 
source, and MOS technology has proved to be suitable for the economical, single chip 
implementation of analog/digital systems. This trend towards mixed analog/digital 
chips, not only for the direct interfacing to the physical world, but also for aiding 
the digital systems to increase their performance will continue in future. The recent 
examples which are the manifestation of this trend are the disk-drive chip[l] and the 
analog assisted microprocessor chip[2]. It is predicted that in few years most chips will 
contain at least some analog circuits in them. 

In the era of VLSI, the analog and mixed-signal designs are becoming popular 


2 



because of the reasons described above. Also the subthreshold region is getting more 
importance day by day. Previously, this region was considered mainly for the calcula- 
tion of the off-stage current in digital ICs because this region was not exploited in the 
application point of view in digital circuits. Now with the advent of analog ICs, there 
are quite a few applications where the subthreshold region of operation is exploited, 
e.g., in digital wrist watches, pacemakers, space electronics, etc. 

The computer simulation of the circuits is one of the crucial factors contribut- 
ing to the success rate of ICs. The circuit simulation is done after the circuit design 
and prior to the fabrication of the chip. For accurate simulation of the electronic cir- 
cuits, the device has to be modeled accurately. The device models developed earlier 
have been proved to be very accurate for long channel MOSFETs. However, for short 
channel MOSFETs, these models are proved to be quite inaccurate. This is due to 
the various short channel effects present in miniaturized devices, which are ignored in 
the long channel MOSFET model development due to their negligible effect on them. 
However, these effects are the dominant ones for small geometry devices. Hence, the 
short channel MOSFET modeling in its subthreshold region of operation has become 
a challenging and important area. 


1.1 The Device Modeling 

The device model is a mathematical representation of the various physical phenomena 
occurring in the device. The device model consists of a set of mathematical equations 
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which are used to predict the behavior of the device accurately under the known excita- 
tions of the various input voltages. The device model equations have a set of parameters 
termed as the model parameters. The model parameters which are used to represent 
the device phenomenon analytically are referred to as the physical parameters, whereas 
the ones which are used in the model equations such that these equations are the nu- 
merical approximations of the various device phenomena are known as the empirical 
parameters. The development of the models whose equations are totally based on the 
analytical representations of the physical phenomena is a difficult task. These mod- 
els are obviously highly accurate, however, they use up too much of CPU time, and 
hence, are inefficient for circuit simulation. On the other hand, the totally empirical 
approach of device model development generating the model having too many empiri- 
cal parameters is computationally less expensive but also is less accurate and does not 
always give a true physical picture of the events taking place within the device. Also, 
the values of these empirical parameters are dependent on the device dimensions and 
specifications such as effective channel length Leff, channel width W, substrate doping 
Nsub, and oxide thickness Tqx- Hence, there is a clear trade-off between the accuracy 
of the model and its complexity. 

The behavior of the circuit is also predicted using the device model, and hence, 
the accuracy in the prediction of the circuit behavior depends upon the accuracy of 
the device model used in the circuit simulation. The circuit simulation (i.e., the pre- 
diction of the behavior of a circuit) may not be adequately accurate even if the device 
model used in the circuit simulation is highly accurate under the DC condition. The 
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various important requirements of the MOSFET model for incorporating it in a circuit 
simulator are discussed in the next section. 

1.2 The Desirable Characteristics of a MOSFET 
Model Suitable for Analog Circuit Simulation 

The device model is used for the bias point calculation in the circuit simulation as 
well as for calculating the small signal quantities such as the transconductance g-m 
(i.e., the ratio of the change in the drain current with respect to the change in the 
applied gate voltage), and the drain conductance Qd (i.e., the ratio of the change in 
the drain current with respect to the change in the applied drain voltage). All these 
small signal quantities along with the drain current should be continuous with respect 
to the terminal voltages. Several benchmark tests have been proposed by Tsividis and 
Suyama[3] for testing the performance of MOSFET models which can be used in analog 
circuit simulation. The characteristics which are to be exhibited by a MOSFET model 
in order to make it suitable for analog circuit simulation are briefly described below. 

1. The model should account for the Drain Induced Barrier Lowering (DIBL) effect 
accurately. This effect is the lowering of the source-channel potential barrier 
in the subthreshold region due to the applied drain-to-source voltage in devices 
having very short channel lengths. This results in a considerable increase in the 
drain cmrent in the subthreshold region of operation. It has neghgible influence 
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in long channel MOSFETs, but ignoring this effect in the model development for 
short channel MOSFETs results in a large error in the subthreshold region drain 
current calculations. 

2. In short channel devices, an increase in the drain voltage results in an increase in 
the subthreshold slope, where the subthreshold slope is the ratio of the change 
in the applied gate voltage to the change in the channel surface potential at the 
midpoint of the subthreshold region of operation, and hence, this affects the value 
of the drain current too in the subthreshold region. This is also a short channel 
effect which has to be considered in the model development in the subthr^hold 
region. 

3. The transconductance Qm becomes a constant in the strong inversion region (as- 
suming that the Id-Vgs characteristics can be expanded into a the first degree 
polynomial). On the other hand, in the weak inversion region, the drain current 
is dependent exponentially upon the applied gate voltage, and, hence, the ratio 
of the transconductance to the drain current (i.e., gm/Io) becomes a constant 
in the weak inversion region. This ratio, which is thus a constant in the weak 
inversion region, decreases as the device moves into the strong inversion region. 
This function (i.e., gm/Io) should be continuous with respect to Vgs- In many 
existing MOSFET models, this ratio becomes discontinuous in the moderate in- 
version region, i.e., at the boundary between the weak and the strong inversion 
regions. This is because the Id-Vcs characteristic is neither exponential nor can 
be expressed by a first degree polynomial in the moderate inversion region[4]. 
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The calculation of the transconductance in the moderate inversion region using 
either the strong or the weak inversion region model results in a large error. With 
the advent of low supply voltage designs, the moderate inversion region has be- 
come extremely important. The discontinuity in the transconductance results in 
a large error in the calculation of the gain of analog circuits. Thus, the continuity 
in Qm should be ensured in order to make the model suitable and accmate for 
analog circuit simulation. 

4. The drain conductance should be continuous with respect to the applied drain 
voltage VjDs- Most of the available models show an abrupt behavior of the drain 
conductance at the boundary between the linear and the saturation region. This 
abrupt behavior of Qd results in a large error in the gain calculation for analog 
circuits. 

The continuities in the g-m/Io ratio with respect to the applied gate voltage and that 
in gd with respect to the applied drain voltage not only give accurate results in analog 
circuit simulation but also malce the convergence in the simulation faster. 


1.3 Literature Review 

A brief overview of some recent MOSFET models is presented in this section. Also, 
the shortcomings of these models are highlighted. 

The first MOSFET model, in which the DIBL effect was modeled by intro- 
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ducing a DIBL parameter S was proposed by 'B:outman[5]. This DIBL parameter was 
purely a fitting parameter, and no analytical expression was given for it. Its value had 
to be found from the measurement data. It had been observed that 6 was a ftmction 
of the oxide thickness, the substrate doping, and the effective channel length. Thus, 
the DIBL effect was not adequately handled by this model. 

The Berkeley Short Channel IGFET Model (BSIM)[6] is one of the most pop- 
ular MOSFET models in recent times. This model is highly accurate, and it is partly 
physical and partly empirical, requiring quite a few fitting paxameters. The subthresh- 
old slope T} in BSIM is given by the following equation[6] 

7? = J?o + (1-1) 

where rjo, r]B, and riD are the model parameters, and the values of these parameters 
have to be found from the measurement data. No physical basis for t)o, tjb, and rjB 
are given. In the drain current model of BSIM, the strong inversion component has 
been matched with the weak inversion component at the boundary between the strong 
inversion and the weak inversion regions in order to restore the continuity in the drain 
current and its first derivative with respect to the applied gate voltage. 

The Unified Charge Control Model (UCCM)[7] is a good model for analog 
circuit CAD (Computer Aided Design). It preserves the continuity in the derivatives 
of the drain current with respect to the terminal voltages at the various operating 
region boundaries. However, its major drawback is that it fails to predict the device 
behavior accurately in the deep subthreshold region, since it does not adequately model 
the effect of the applied drain voltage on the subthreshold slope for very short channel 
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The Physically based Continuous short-channel IGFET Model (PCIM)[8] is 
another analytical model for short channel devices. This model is described by a single 
drain current equation, valid for both the weak inversion and the strong inversion 
regions. This model maintains the continuity in the drain conductance as well as in 
the capacitances with respect to the terminal voltages at the various operating region 
boundaries. The DIBL effect has been modeled in this model. Since this model is based 
on the device physics, it automatically accounts for the changes in the behavior of the 
device with the changes in the device parameters such as the oxide thickness, substrate 
doping, temperature, and channel length. However it fails to model the behavior of 
the device in the deep subthreshold region. 

A modified subthreshold slope model has been proposed by Kang et al.[9]. In 
this model, the values of the model parameters are determined from the threshold volt- 
age data. This model adequately handles the DIBL effect by successfully expressing 
the subthreshold slope in terms of the threshold voltage parameters, after account- 
ing for both the body and the DIBL effects. However, this model does not properly 
describe the effect of increase in the subthreshold slope with increase in the applied 
drain voltage. In order to develop a new and more accurate subthreshold slope model, 
this model[9] can be modified to make the subthreshold slope a function of the applied 
drain voltage. As mentioned before, the subthreshold slope equation in BSIM is purely 
empirical, where the values of the fitting parameters strongly depend upon the device 
specifications, such as Tox, Ngub-, Lg//, etc. This subthreshold slope model in BSIM 
can be replaced by the new subthreshold slope model developed in this work after 
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can be replaced by the new subthreshold slope model developed in this work after 
modifying the subthreshold slope model proposed by Kang et al.[9] in order to make 
BSIM more accurate. ■ 


1.4 A Brief Description of This Work 

In this work, a new analytical model for the subthreshold slope has been developed. 
This model is strongly based on the device physics, and the subthreshold slope depen- 
dencies on the effective channel length, the body voltage Vsb and the applied drain 
voltage Vds have been incorporated analytically into it along with the DIBL effect. 

The subthreshold slope model proposed by Kang et al.[9] has been modified 
to come up with a new subthreshold slope model. The original subthreshold slope 
model [9] has been developed from the capacitor equivalent circuit of short channel 
MOSFETs in the subthreshold region of operation. In [9], the values of the capacitors 
in the circuit are analytically expressed in terms of the channel length after assuming 
that the parameters appearing in the expressions for the capacitances are independent 
of the drain voltage. Therefore the subthreshold slope in that model[9] becomes in- 
dependent of the drain voltage. However, from the experimentally observed results, 
it has been observed that the subthreshold slope 77 is a strong function of the drain 
voltage, and it increases considerably with an increase in the drain voltage particularly 
for small channel length devices. Basically, in the subthreshold slope model proposed 
by Kang et al.[9], it was assumed that in the subthreshold region, the drain voltage 
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does not change the physical picture of the MOSFET channel depicted in terms of 
the capacitor equivalent circuit. However, it has been observed in this work that the 
point of the channel where the suff^e potential is minimum is a strong function of 
the drain voltage, the body voltage, and the effective channel length. This implies 
that, in the subthreshold region, the drain voltage, the body voltage, and the effective 
channel length modulate the values of the capacitors of the equivalent circuit given by 
Kang et al.[9]. This makes the subthreshold slope a function of these three parameters. 
We have analytically modeled the effects of these parameters on the minimum surface 
potential point, in order to develop an improved analytical subthreshold slope model 
which accounts for the dependencies of the subthreshold slope on the effective channel 
length, the drain voltage, and the body voltage. 

In order to develop a drain current model, the subthreshold slope model and the 
threshold voltage model used in the drain current model developed by Deshpande and 
Dutta[10] have been replaced by the new subthreshold slope model developed in this 
work and the BSIM threshold voltage model respectively. The results of the improved 
subthreshold slope model and modified drain current model developed in this work 
have been compared with the experimentally obtained results reported in literature for 
MOSFETs having effective channel lengths ranging from 1.63 /rm to 0.075 /xm. The 
simulated results show a good match with the experimentally obtained ones. 
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1.5 The Organization of This Report 


The development of the improved subthreshold slope model based on the BSIM thresh- 
old voltage model is discussed in detail in Chapter 2. This chapter also includes the 
development of the modified drain current model by replacing the subthreshold slope 
model and the threshold voltage model used in the drain current model developed 
by Deshpande and Dutta[10] by the improved subthreshold slope model developed in 
this work and the BSIM threshold voltage model respectively. The extraction proce- 
dure for the values of the various parameters used in the improved subthreshold slope 
model from the threshold voltage data is presented in Chapter 3. Also the simulated 
results of the improved subthreshold slope model and the drain current model have 
been presented in this chapter. The results from the experimentally obtained data 
for various short channel MOSFETs have been given along with the simulated results 
obtained from our model for the sake of comparison. The results show an excellent 
match. The summary of the entire work has been presented in Chapter 4 along with 
the conclusions. Also, the shortcomings of this model along with the future scope for 
improvement have been discussed in this chapter. 
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Chapter 2 


Model Development 


2.1 The MOSFET Structure 

The basic three dimensional structure of an n-channel MOSFET is shown in Fig.2.1, 
along with the axes chosen. The channel length is along the x-direction, and the 
channel width is along the z-direction. The two regions (typically few tenths of a 
/im deep), which axe diffused in the p-type substrate (or body) are called the source 
and the drain; and the inversion layer, which is formed in between these two regions is 
generally referred to as the channel. The structure also has a very thin layer of oxide 
(Si02) grown on the p-type substrate. The oxide thickness is denoted by Tqx-, and it is 
typically several hundreds of A. The oxide is followed by a gate material (shown by the 
shaded region in the figure), which can be either metal or polysilicon'. The drawn gate 
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length L is the length of the oxide which separates the ri^ drain and source diffusions. 
This is also known as the metallurgical channel length. When the drain and the source 
are actually diffused into the substrate, the effective channel length Lg// becomes less 
than the drawn length due to the lateral diffusion of the drain and the source dopants. 
The dashed line in the figure is the depletion region, whose width along the y-direction 
is referred to as Xdej>- Obviously, this depletion width is a function of position along 
the length of the channel. 



Body 


Figure 2.1: The Structure of an n-MOSFET. 
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2.2 The Subthreshold Slope Models 


Theoretically, the subthreshold slope, represents the changes in the surface potential xpa 
with respect to the applied gate voltage Vqs in the subthreshold region of operation. 
Also, the threshold voltage depends upon the surface potential variation with respect to 
the gate voltage. The changes in the geometry of the MOSFET, the substrate doping 
Nsub, and the oxide thickness Tqx affect the surface potential variation with respect 
to the gate voltage. This, in turn, affects the value of the subthreshold slope and the 
threshold voltage. Therefore, accurate subthreshold slope modeling is closely related 
to precise threshold voltage modeling. Actually, the subthreshold slope model can be 
developed from an accurate threshold voltage model. 


2.2.1 The Classical Subthreshold Slope Model 


A brief introduction to the classical subthreshold slope model[9] is given here. The 
various parameters used in this model are obtained from the threshold voltage data. 
The basic definition of the subthreshold slope or the ideality factor rj is given by [11] 


V 


5Vgs 


( 2 . 1 ) 


where Vgs is the gate-to-source voltage (henceforth referred to as the gate voltage), 
and tps is the surface potential. 


The ideality factor for long channel MOSFETs can be represented by[5] 


T] = 1 + 


Cb 

Cox' 


( 2 . 2 ) 
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where is the depletion capacitance per unit area, and Cqx is the oxide capacitance 
per unit area. 

It has been found that for long channel devices, the results obtained from this 
model match weU with the measurement results [5]. However, this model cannot be 
used for short cheumel devices due to various second-order effects. For such devices, 
the ideality factor model has to properly consider the various short channel effects. 


For short channel MOSFETs, the threshold voltage Vts is given by [6] 

Vts = VpjB -H 24>f + + ^SB — r(2^jr + Vsb) — ctVds, (2-3) 


where Vps is the flat band voltage, 4>f is the bulk potential, Vsb is the applied source- 
to-body voltage (henceforth referred to as the body voltage), Vbs is the applied drain- 
to-source voltage (henceforth referred to as the drain voltage), 7 is the body effect 
coefficient, T is the nonuniform doping effect parameter, and a is the Drain-Induced 
Barrier Lowering (DIBL) parameter. 


Since the subthreshold slope 77 is defined as (^fvcs) Eqns.(2.2) and 

(2.3), the expression for can be obtained by substituting Vqs for Vts, tps for 2^5-, 
and taking the derivative of Vqs with respect to t/'j. Thus, can classically be 
expressed as [11] 


Cb ^ 7 

Cox 2^/l.b(j)F + 


(2.4) 


Thus, from Eqns.{2.2) and (2.4), the classical model of the subthreshold slope ij can 
be given by 
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(2.5) 


7 

^ 2\/l.5(^ir + VsB 

It should be noted that the DIBL effect, which is extremely important for short channel 
MOSFETs, has been neglected in this model [Eqn.(2.5)]. In this work, an analytical 
subthreshold slope model is developed taking into account the DIBL effect in a com- 
prehensive manner. 


2.2.2 The DIBL Effect 


In the subthreshold region, there is a capacitive coupling present between the 
channel and the source/drain, which can be represented by two coupling capacitances 
Csc and C dc respectively. These capacitances are in addition to the already existing 
capacitances Cox and Cb, and are given by 


Csc = 



( 2 . 6 ) 


and 


^DC r 1 


(2.7) 


where Csi is the permittivity of Si, k is the proportionality constant, Leff is the effective 
channel length, and Xo is the distance along the channel (measured from the source) 
where the surface potential is minimum. Thus, the equivalent circuit of a short channel 
MOSFET at subthreshold can be given as shown in Fig.2.2. 

Therefore, the surface potential changes not only due to the gate voltage Vgs 
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but also due to the drain voltage Vds- The parameter o appearing in Eqn.(2.3) models 
this effect. The values of these coupling capacitances {Cdc and Csc) increase as the 
channel length is reduced. The value of the DIBL parameter a also increases as the 
channel length is reduced, which implies that the DIBL effect is more pronounced in 
short channel MOSFETs. This makes it important to model this effect accurately with 
the help of the parameter a. This parameter is a function of the effective channel 

Gate 



Source 



C 


ox 


Drain 



Substrate 

Figure 2.2: Capacitor equivalent circuit of a short channel MOSFET at subthreshold. 

length Leff and the body voltage Vsb- In order to find the variation in the value of a 
with respect to the term ■\/ 24 >f + ^sb for a particular channel length, AVts(AVds (i-e-, 
the change in the threshold voltage Vts with respect to the change in the drain voltage 
Vos) can been plotted with respect to \/2((>f + Vsb for that particular channel length. 
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Figure 2.3 shows the plots of {AVts/ AV^ s) / y/24>F + Vsb ^ a function of y/24>F + Vsb 
for four different values of the effective channel lengths (i.e., for Leff = 0.33 /im, 0.43 
Ijm, 0.63 jum, and 1.43 yum). These data are obtained from the results reported by 
Lee and Min[13]. From this plot, the values of a can be obtained for various channel 
lengths. These results are described in the next chapter. 



/ 172 

V2(t)p +VsB (Volts) 

Figure 2.3; {AVts/ AV^ s) /{2(pF + VsbY^^ versus {2(I)f + VsbY^"^ for different channel 
lengths. 

The DIBL parameter a has been expressed as [9] 

a = s^2<Pf + VsB — S{2<j)F + Vsb), (2-8) 

where s and S are the model parameters, which are used to model the behavior of a 
with respect to Vsb- Thus, one can obtain the values of 7 , F, s and S by the method 
of curve fitting on the threshold voltage data. The results obtained from this type of 


19 



parameter extraction are described in the next chapter. 


The channel potential can be calculated by the capacitance-coupling theory 


as [9] 


Vs 


WsO[^GO^ySB,Vso,yDo)+ rv , ^ ^ 

^‘OX + "T Odc + L^SC 


(2.9) 


In Eqn.(2.9), ipg has been expressed by expanding the Taylor series around ipso 
— l.S^F + VsB (which is the approximate middle point of the subthreshold region), 
and keeping only the first-order terms. Since the corrections in the values of Vts and 
due to the capacitances Cue and Csc are quite small, it is enough to keep only the 
first-order terms in the Taylor series expansion. 


2.2.3 The Subthreshold Slope Considering the DIBL Effect 

The ideality factor 77 is obtained by partial differentiation of Eqn.(2.9) with respect to 
Vg, and it can be expressed as [9] 


77 = 1 - 1 - 


Cb + Csc + Cdc 


a 


ox 


or 


77 1 -f 


Cb 

Cox 


+ 


2Cdc 

Cox 


( 2 . 10 ) 
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In obtaining Eqn.(2.10), the subthreshold slope 77 is evaluated at Vbs = 0, and, there- 
fore. the capacitances Csc and Cdc have the same value, since the minimum surface 
potential point is located at the middle of the channel, i.e., at Xo = Lsffl2. Thus, from 
Eqns.(2.6) and (2.7), one gets 


Cdc{Vds = 0 ) = Csc{Vds = 0 ) = 2esik/Leff- ( 2 - 11 ) 


By partial differentiation of Eqn.(2.9) with respect to Vjos, and neglecting Cdc 
and Csc terms appearing in the denominator, the following relationship can be 
obtained [9] 

a — 

Cb + Cox + Cdc + Csc ' 


or 


(7 


Cdc 

Cb + Cox 


( 2 . 12 ) 


The term Cdc appearing in Eqn.(2.12) is the drain-channel coupling capaci- 
tance for large values of Vds [9]. From Eqn.(2.7), it can be clearly seen that Cdc is 
inversely proportional to Le// — Xo- Though the minimum surface potential point Xo 
is known to be a complex function of the channel length, the oxide thickness, and the 
drain voltage [14], it has been stated that for large values of Vds-, if becomes indepen- 
dent of Vds and converges rapidly to a value equal to Leff/^ [15]. Therefore, for large 
values of Vds, Cdc can be expressed as [9] 

CDc{largeVDs) = l'f'^- (2.13) 

2Leff 

Thus, from Eqns.(2.11) and (2.13), one can write 

= = (2.14) 
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Similarly, from Eqns.(2.12) and (2.14), one gets 


Cdc{Vds = 0 ) = 


4cr(CB + Cox) 
3 


(2.15) 


Thus, from Eqns.(2.4), (2.10), and (2.15), one gets the final expression for the 
ideality factor, which is given by [9] 


In this derivation, it has been assumed that the minimum surface potential 
point lies at a distance equal to Leff/S from the source. This assumption makes 
the ideality factor independent of Vos, which is incorrect, as has been proved by the 
experimental data on short channel length devices reported in literature, which show 
changes in the subthreshold slope rj with a change in the applied drain voltage Vos- 


The development of the improved subthreshold slope model considering the 
effect of the applied drain voltage Vds is described in the following section. 


2.2.4 The Improved Subthreshold Slope Model 

This section describes the simulation algorithm, which accounts for the dependence 
of r/ on Vos- Iii reality, the minimum surface potential point Xo is a function of Vds 
and Leff- Thus, if these variations can be modeled in the expression for Xo-, then 
the expressions for Cdc aiid Csc would also be altered, with each having a specific 
dependence on Vds and Leff- And since t? is a function of Cdc and Csc-> if will have 
an explicit dependence on Vds and Te//, which previous authors have failed to model. 
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The algorithm is as follows. 


1. Evaluate the surface potential at each point of the channel for various values of 
Leff and Vds- 

2. Find the values of Xo (i.e., the distance of the minimum surface potential point 
from the source) for different values of Lgfj and Vos- 

3. Use these values of Xo (having strong dependence on Vds and Tg//) for finding 
the values of Gdc and Csc- 

4. Then, compute 77 as a function of Vds and Lg/f from the expression given earlier. 


The details are given below. 

The built-in voltage of the source-substrate junction is given by 

Hi = (2.17) 

where Vj is the thermal voltage (= kT/q), Ngub and Nd are the substrate and the 
source doping concentrations respectively, and n* is the intrinsic carrier concentration. 
The depletion layer thickness at the source end (at threshold) is given by [11] 

Xdep = + VsB)/qNsub- (2-18) 


The surface potential '^^(x) as a function of Vds and Lgff (at Vsb = 0) can be given 
by [14] 

, / ^ , /X. rr , .sinh(xfl) , .sinhUL — x)/l] 

'ips{^) = ipsL + {Vbi + VDs-AL)^.^^f^j^^j^^+{Vbi-'ilJsL) ■’ ( 2 - 19 ) 
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where 'ipsL is the long channel surface potential, given by [14] 

‘ipsL = Vgs — Km + 2^i?, (2.20) 

with VthQ being the long channel threshold voltage, given by[14] 

Vm = Vfb + 20 ^, ( 2 . 21 ) 

^ox 


where Tqx is the oxide thickness, and cqx is the permittivity of the oxide. The 
characteristics length 1 is defined by [14] 



.The depletion layer thickness Xdep at threshold has been assumed to be uniform along 
the channel. However, in reality, it is a function of V^s and Leff, particularly so 
for short channel length devices. Therefore, in order to compensate for this invalid 
assumption (of constant Xdep), another factor rjf has been introduced in the expression 
for I [Eqn.(2.22)], which is basically a fitting parameter called as the fudge factor. 
Xdep/Vf inay be taken as an average of the depletion layer thickness along the channel. 
The fudge factor rjf itself may be a function of V^s, however this is a second order 
effect, and hence it is neglected [14]. 


It has been observed that for long channel devices in subthreshold, the sur- 
face potential is almost constant throughout the channel. Hence, the value of rjf can 
be taken to be equal to unity. However, for short channel devices (having submicron 
channel lengths), the portion of the channel where the surface potential stays con- 
stant is quite small, and, hence, the value of rjj is correspondingly higher. From the 
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And finally, the improved expression for the ideahty factor developed in this work can 


be given by 

b = T/ - r) + 4(1 - )cr) . (2.26) 

V 2-v/l.O0f’ + VsB ) V -^e// / 

Equation (2.26) is the improved subthreshold slope model developed in this work, and 
will be used for estimating the subthreshold region drain current, as illustrated in the 
next section. The results are presented in the next chapter. 


2.3 The Subthreshold Region Drain Current Model 


The subthreshold region drain current expression for a short channel MOSFET is given 
by [10] 


Leff 2'0° 


(2.27) 


where iin is the electron mobility, ip^ is the short channel surface potential 
given by [10] 

= Vgs — Vtho-^‘^4’F + VsB, (2.28) 

Vtho is the long channel threshold voltage [Eqn.(2.21)], Vts is the short channel thresh- 
old voltage [Eqn.(2.3)], p is the reciprocal of the thermal voltage Vp, and r) is the 
subthreshold slope [Eqn(2.26)]. 
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2.4 The Strong Inversion Region Drain Current 


Model 


The strong inversion region drain current model is given by Sheu et al.[6]. This model 
has been modified by Deshpande and Dutta[10]. The modification algorithm of Desh- 
pande and Dutta[10] is briefly reproduced here. The strong inversion component of 
the drain current Inn in the linear region is given by [6] 


fyarVBs - ivis) 

Leff 1 + OVds + UqVqt (1 + ^^^^-Vds) 


{Vgt > 0, 0 < VjDs < Vdsat) 

(2.29) 


where [4] 


Vqt — Vgs — Vts, 


(2.30) 


a 


1 + 


91 

2y/2(j)F + VsB ’ 


(2.31) 


and 


1.744 + 0.8364(2<^ir + Vsb) ’ 


(2.32) 


Uo is the vertical field mobility degradation coefficient, Ui is the velocity saturation 
coefficient, and 6 accounts for the source and drain parasitic series resistances [6]. 

The parameters 7 and T model the non-uniform doping profile [4]. In addi- 
tion to modeling the DIBL effect, the parameter a accounts for the channel length 
modulation eflFect as well[6]. 
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The drain current expression in the saturation region {Vgs > and Vos > 
Vdsat) is given by [6] 


I ^ C 

Leff 1 + OVds + UqVqT ‘2-Oi K 


(2.33) 


where 


^ _ 1 + Uc + \/\ + 2Uc 
2 


(2.34) 


with 


U\ Vqt 
Wc= r 

.be// OL 

The drain-to-source saturation voltage Vdsat is given by [6] 


(2.35) 


V 


Vgt 


DSAT — , 


aVK' 


(2.36) 


An additional parameter li had been introduced by Deshpande and Dutta[10] 
in order to take into account the channel length modulation effect in the saturation 
region drain current expression, with the final expression given by [10] 


T' ^ W C 

Lef / ""^ 1 + eVns + UoVgt 2a K{1 - h (Vds - Vdsat) ) ' 


(2.37) 


However, the DIBL parameter a partially accounts for the channel length modulation 
effect [6]. Therefore, the threshold voltage model [Eqn.(2.3)] used here already accounts 
for this effect. Hence, one need not incorporate the parameter li separately for it. 
Therefore, Eqns.(2.27), (2.29), and (2.33) give the expressions for the drain currents 
for subthreshold, linear and saturation regions respectively. 
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2.5 Approaches for Joining the Various Regions 


of Operations of the MOSFET 

The discontinuities in the drain current (Id), the drain conductance (g^), and the 
transconductance {gm) with respect to the terminal voltages at the boundaries between 
the various regions of operation create large errors in circuit simulation, especially for 
analog circuits [3]. The most basic requirement for a MOSFET model fit to be used 
in analog circuit simulation work is the continuity of Id versus Vos characteristics. 
The discontinuities occur at two points in the operating range, one at the boundary 
between the weak and the strong inversion regions (i.e., at Vgs = Vts), and the other 
at the boundary between the linear and the saturation regions (i.e., at Vds = Vdsat)- 

The discontinuity of /£> with respect to Vqs at the onset of strong inversion 
is due to a total disappearance of the strong inversion component of /x? in the weak 
inversion region (i.e., for Vos < Vts). The improper models used for the weak and the 
strong inversion regions can also lead to a discontinuity in Id with respect to at 
the boundary between the weak and the strong inversion regions [3]. The discontinuity 
in the derivative of In with respect to Vds occurs at the boundary between the linear 
and the saturation regions. In the saturation region {Vos > Vts, and Vds > Vdsat), 
Id is almost constant with respect to Vds (neglecting the channel length modulation 
effect) with a value almost equal to that obtained for Vds = Vdsat, whereas in the 
linear region (Vos > Vts, and Vds < Vdsat), Id varies almost linearly with Vds- This 
abrupt change of characteristic for Id from the linear to the saturation region makes 
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the derivative oi Id discontinuous with respect to Vds- 


2.6 Continuity in the Drain Conductance with Re- 
spect to the Drain Voltage 


As mentioned before, the discontinuity in the drain conductance {gd) with respect 
to the applied drain voltage ^£>5 arises at the boundary between the linear and the 
saturation regions. The drain current is a strong function of the drain voltage Vds 
in the linear region, whereas in the saturation region, it becomes essentially constant. 
This implies that the drain current Id has completely different dependencies on the 
drain voltage Vds in these two regions. This change in the behavior oi Id with respect 
to Vds occurs at the boundary between the linear and the saturation regions, and it is 
abrupt. Therefore, Vds can be viewed as if it saturates to Vdsat in an abrupt manner 
at the onset of saturation (i.e., at Vds = Vdsat)- In order to remove this discontinuity, 
Deshpande and Dutta[10] have made Vds approach Vdsat gradually. This has been 
accomplished by replacing Vds in Eqn.(2.29) by “effective drain-to-source voltage Vdax' ■ 
The expression for Vdsx is given by [10] 


Vdsx = Vdsat 1 + 


f ln[l + exp(A{l 


_ yps 
Vdsat 


))]^ 


ln[l + exp{A)\ j 


(2.38) 


where the coefficient A is a fitting parameter, and the sharpness in the transition of 
Vdsx to Vdsat depends on its value. 
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Equation (2.38) has been simulated and plotted in Fig.2.4, which clearly shows 
that the value of V^sx gradually approaches Vdsat beyond the point Vds = Vdsat- In 
the linear region, the value of Vdsx is almost equal to the value of Vos- The behavior 
of Vdsx around the transition point (i.e., at Vos = Vos at) depends upon the value of 
A. From Fig.2.4, the effects of the value of A on the transition of the curve can be 
observed. For curves having sharp transitions, the values of A are high, whereas for 
gradually changing curves, the values of A are low. It has been observed that for very 
short channel devices, the value of A is low because of the widening in the transition 
width [10]. The width of the transition from the linear to the saturation region is low 
for long channel devices, therefore the value of A is correspondingly higher. 



Figure 2.4: The normalized effective drain-to-source voltage {V^sx/Vdsat) versus the 
normalized applied drain-to-source voltage (Vos/Vdsat) for various values of A. 
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From Eqns.(2.29) and (2.38), a single expression has been derived by Desh- 
pande and Dutta [10] for the strong inversion region drain current, which can be given 


by 


Leff °^(l+mS + UoVGT) 


(2.39) 


Equation (2.39) maintains the continuity of the drain conductance with respect 
to the applied drain voltage Yds at the boundary between the linear and the saturation 
region. These results are shown in the next chapter. 


2.7 Continuity in the Transconductance (g^) with 
Respect to the Gate Voltage 


The transconductance (ffrn) should be continuous with respect to Vos at various op- 
erating region boundaries. Discontinuity in g^n occurs at the boundary between the 
weak and the strong inversion regions (i.e., at Vos = Vts)- This discontinuity is re- 
moved by Deshpande and Dutta [10] by making the weak inversion component of the 
drain current slowly saturate to some value as Vas is increased beyond Vxs- This is 
accomplished by modifying the weak inversion component of the drain current as [10] 


^ sub 


I sub 

I sub ” 1 “ I limit 


where the term I sub is given by [10] 


(2.40) 


•-suh 


n V'‘ f, _ .-tlVos)MVss-VT,)h 


( 2 , 41 ) 
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where M is a fitting parameter, which is used to describe the moderate inversion re- 
gion. and the value of this parameter is empirically chosen to achieve best fits in the 
subthreshold region characteristics with minimum effect on the strong-inversion region 
characteristics [6]. lumit is given by [10] 

~ ^sub\vQx—^' (2.42) 

The higher order continuity in the transconductance with respect to the applied gate 
voltage Vqs is restored by Deshpande and Dutta [10] by making the strong inversion 
component I^tr gradually decrease to zero as Vqs decreases below Vts- This gradually 
decreasing drain current component Igtr in the weak inversion region should be much 
smaller than the actual weak inversion component. In order to accomplish this, the 
voltage Vqt is modified to an effective voltage Vgtx such that [10] 

Vgtx = Si Zn ^1 -1- ^ ^ , (2.43) 

where Si is a fitting parameter. This equation has been simulated and the normal- 
ized values of Vgt^ have been plotted in Fig.2.5 for two different values of the fitting 
parameter Si (O.SV^ and Vr) as a function of the normalized values of Vct- In the 
strong inversion region {Var > 0); Vgtx approaches Vgt- On the other hand, in the 
weak inversion region {Vgt < 0), the values of Vgtx are much less than Vgt- Thus, 
the behaviour of the Vgtx in the Fig.2.5 implies that the strong inversion component of 
the drain current in the weak inversion region holds some nonzero value which is much 
much less than the weak inversion component value. 
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Figure 2.5; Normalized values of Vgtx for two different values of Si {0.5Vr and Vr). 


2.8 The Complete Model 


The total drain current is the sum of the weak and the strong inversion components, 
which can be given by 

Id — I sir + I sub ) (2.44) 


where Istr is the combined strong inversion component of the drain current, valid for 
both the linear and the saturation regions, given by [10] 


^str 


W IXn VgtxV^isx - 

Leff °^{l + eVDS + UoVgtx) (1 + /) Vd,,) ’ 

e/ j 


(2.45) 


and the weak inversion component I sub of the drain current is given by Eqn.(2.40). 
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The fitting parameters k' , Si, d, Uq, Uu A, 7, T, a, rjfo are to be extracted 
from the measurement data. The parameters cr, 7 and F can be extracted from the 
threshold voltage data. The values of the parameters 77/0 and k' are to be extracted 
from the measurement data of the subthreshold slope. The results simulated from this 
model are given in the next chapter. 



Chap;;er 3 


Model Results 


In this chapter, the simulated results of the models described in the previous chapter are 
presented. Initially, the simulated results for the subthreshold slope model described by 
Kang et al.[9] are presented, and the inherent flaws in that model are highlighted. Next, 
the simulated results of the improved subthreshold slope model developed in this work 
are presented. The parameters required for this model are extracted from the threshold 
voltage data. The current-voltage characteristics simulated from this improved model 
are also given. The model results are compared with the experimental data reported 
in literature for devices having channel lengths less than hundred nanometers, and the 
results show a good match. 
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3.1 Simulation Results of the Subthreshold Slope 
Models 

In this section, the simulated results of the subthreshold slope models described in 
the previous chapter are presented. For this simulation, MOSFETs having effective 
channel lengths of 1.43 yum, 0.63 yum, 0.43 yum, and 0.33 yum are considered. The oxide 
thickness is 100 A and the substrate doping is 1.6x10^’’ per cm^. We have the Id versus 
Vgs characteristics for these devices for Vsb = 0, 1, 2, 3, and 4 V, and for Vds = 0.1 
V and 3.0 V[13]. These characteristics are plotted in Figs.3.1(a) and (b) for channel 
lengths of 0.33 ^m and 0.63 yum respectively, for different values of the body voltage 

3.1.1 The Subthreshold Slope Model based on the Threshold 
Voltage Data (without considering the effect of Vds on 

n) 

In this section, the method of extraction of the values of the parameters s, S, 7, 
and r from the threshold voltage data for the four devices having different channel 
lengths are presented. From these parameters, the values for the subthreshold slope 
are determined from Eqn.(2.16). This model for the subthreshold slope is independent 
of Vdsi as explained in the previous chapter. 

In order to obtain the values of the parameters s and S, the term 
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(a) 



(b) 

Figure 3.1: Experimentally obtained Id-Vgs characteristics of the MOSFETs having 
effective channel lengths of (a) Leff = 0.33 /xm and (b) Lg// = 0.63yLtm.[13]. 
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{/\Vts/ ^Vds)/ \/24)f + VsB has been plotted as a function of ^/2(j)F + Vsb in Fig.2.3. 
The term AV^s is defined as [6] 


APts = Vts{Vds = 0.1 V) - Vts{Vds = 3.0 V). (3.1) 


The threshold voltage for V^s = 0.1 V has been obtained by using the following 
rule of thumb [7, 11] 


Vts (Vos = 0.1 V) = Vqs 


{g-m —dmmax/^) 


(3.2) 


where the term gm is the transconductance, defined as the first derivative of In with 
respect to Vos, given by 


9m — 


SI 

SVr 


D 


GS 


(3.3) 


and 9mma:c is the maximum value of gm- The gm versus Vqs characteristics have been 
plotted in Figs. 3. 2(a) and (b). From these figures, it can be observed that proper 
maxima exist for Vbs = 0.1 V, but no such maxima exist for Vds = 3 V. This is due to 
the saturation effect which occurs for Vds = 3 V. In saturation, the transconductance 
gm varies almost linearly with the applied gate voltage Vos- Therefore, due to this 
effect, the method given in Eqn.(3.2) for obtaining Vts can not be applied for Vds = 
3 V. Hence, for these cases, the values of Vts have been obtained from the method of 
constant current threshold voltage, given by [11,12] 


Vts\vds = ^-OV — = (3-4) 

where Id = 5 gA is chosen to be the reference value for the drain current. 

As stated by Kang et al.[9], for relatively longer channel MOSFETs (e.g., for 
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Figure 3.2: Transconductance Qm versus Vos for (a) Le/f — 0.33 /xm and Vds — 0.1 V 
and (b) Le// = 1-43 /xm and V^s = 3.0 V. 



Leff = 1.43 Aim), the term {AVts/AVds)/V‘^^f + Vsb is observed to be almost con- 
stant, independent of ■\/2i(j)p -t- V53. It indicates that is directly propor- 

tional to the factor y/2<j)F + Vsb in such MOSFETs. It can be observed from Eqn.2.3 
that the factor AVts/AVds has the same functional behavior as that of a. The factor 
AVts/AVds has been modeled by the expression given by Eqn.(2.8). The values of 
the parameters s and S have been obtained from curve fitting of Eqn.(2.8) with Fig.2.3. 
It has been observed that the value of s (and ultimately that of a) increases as the 
channel length is reduced. This implies that the drain-induced barrier lowering effect 
becomes more pronounced as the channel length is reduced. The behavior of AVts 
follows the same pattern as suggested by Liu et al.[14]. 

The values of the parameters 7 and F have been found by curve fitting of 
Eqn.(2.3) with the Vps data for V^s = 0.1 V. The extracted values of the model 
parameters 7, F, s, and S are listed in Table 3.1, and they match well with those 
reported by Kang et al.[9]. These extracted values of the model parameters 7, F, s, 
and S are used in Eqn.(2.16) in order to calculate the ideality factor ry as a function of 
Vsb- These characteristics are shown by continuous lines in Figs.3.3(a)-(d) for effective 
channel lengths of 1.43 ^um, 0.63 /im, 0.43 /xm, and 0.33 respectively . 
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Table 3.1 


The extracted values of the parameters used in the subthreshold slope model 
given by Eqn.(2.16) for the MOSFETs having effective channel lengths of 
0.33 jjLva., 0.43 ^m, 0.63 /^m and 1.43 ^m. 


Channel Length (/im) 
1.43 
0.63 
0.43 
0.33 


7 r 

0.550498 -0.027129 
0.590389 -0.005396 
0.604011 0.024737 

0.509722 0.023842 


s (y-V2) s (y-^) 
0.014307 0.004690 

0.029623 0.007481 

0.045488 0.009397 

0.067303 0.015221 


The experimental values of the ideality factor 77 are obtained by the following 
expression 


/J_\ f6ln{lD)^ ^ 

vV^/ ^ SVqs ) 


(3.5) 


Thus, these values are obtained by a first order approximation of the ln(Ji5) versus 
Vqs characteristics for Vds = 0.1 V in the subthreshold region of operation, (i.e., for 
applied Vqs < Vts)- The subthreshold slope (ideality factor) values obtained from the 
measurement data[13] for Vds — OT ^ shown in Figs.3.3(a)-(d) by the diamond 

symbols. The simulated results are observed to be overestimating the ideality factor, 
particularly so for small channel lengths (i.e., for all the channel lengths below 1.43 
/xm). These results are expected since the dependence of 77 on Vsb was not modeled 
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properly by the earlier authors[9]. 


It has been assumed by Liu et al.[14] that the point Xo along the channel 
where the surface potential is minimum converges rapidly to Lg ^^/3 from the source 
end for higher values of Vos- Hence, the ideality factor 77 in its final expression comes 
out to be independent of the drain voltage V£) 5 [ 3 ]. However, this assumption may 
not quite hold true for small channel length devices, since for these cases, the applied 
drain voltage modulates the surface potential profile quite markedly, and, hence, the 
minimum surface potential point may not converge to Le///3, but instead, may become 
a strong function of Vds- Hence, an improved model for 77 is needed which explicitly 
accounts for the effect of V^s- This improved model has been described in detail earlier 
in the previous chapter, and the results simulated from that model are presented here. 

3.1.2 Simulation Results of the Improved Subthreshold Slope 
Model 

In this section, the simulated results of the subthreshold slope model given by Eqn.(2.26) 
are presented. In the previous chapter, it has been shown that the minimum surface 
potential point Xo is a function of the applied drain voltage V^s- And since the sub- 
threshold slope (Eqn.(2.26)) is a function of Xo, hence, this model is implicitly Vds 
dependent. 

The surface potential 'ips along the channel has been calculated with the help 
of Eqn.(2.19). Then, Xdep is calculated using Eqn.(2.18) taking Vsb = 0 V. The 
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Figure 3.3: Results of the subthreshold slope model [Eqn.(2.16)] at Vds - 0.1 V for 
(a) Leff = 1.43 (b) Le// = 0-63 (c) L^jf = 0.43 fim, and (d) Le// = 0.33 /xm. 



characteristic length I is found from this value of using Eqn.( 2 . 22 ). For this 
simulation, the value for the fudge factor 77 / has been taken to be equal to 1 . 0 , which 
may not be exactly true. The value of the long channel threshold voltage Vtho is 
calculated from Eqn.( 2 . 21 ). The calculated values of the channel surface potential 'ipg 
are plotted with respect to the position along the channel in Fig. 3.4 for different channel 
length devices (i.e., for Lgff = 0.33 ^m, 0.43 /j,m, and 1.43 ^m). The intention of this 
activity is to check whether the Xo point really lies at Leff/3 location irrespective of 
Vds and L^ff. From this figure, it can be clearly seen that for relatively longer channel 
MOSFETs (e.g., for Le/f = 1.43 /xm), the surface potential is almost uniform along the 
channel. This implies that the value of Xo is close to Leff/2 instead of Leff/3 which 
was assumed in the derivation of the subthreshold slope [Eqn.(2.16)] by Kang et al.[ 9 ]. 



0 0.2 0.4 0.6 0.8 1 

Normalized position along the channel (x/Lef^ 

Figure 3 . 4 : The variation of the surface potential along the channel for various 
channel lengths at Vos = 3.0 V. 
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The values of Xq/ L^ ff are found for VgB — 0 V from the closed form expression 
given by Eqn.(2.24), and these are plotted as a function of the effective channel length 
Leff in Figs. 3. 5(a) and (b) for Vds = 0.1 V and 3.0 V respectively. From these figures, 
it can clearly be seen that Xo does not converges to particularly so for V^s = 

0 V, which is expected. Also, for higher values of Vr,s (e.g., for Vds = 3 V), the values 
of Xo/Leff are much greater than 0.33 (the assumed value in [9,15]) for channel lengths 
of 0.63 urn and 1.43 fxm. In any case, these values are always higher than 0 . 33 , which 
confirm our earlier prediction of this kind of behavior. 

As explained in the previous chapter, the fudge factor rjf appearing in the 
expression of the characteristics length I [Eqn.( 2 . 22 )] is a function of both Leff and Vsb- 
From a comparison between the results obtained from the subthreshold slope model 
[Eqn.(2.26)] and the measured values of the subthreshold slope from the experimental 
data, the values of the fudge factor rjf have been found for each channel length and 
Vsb- It has been observed that this fudge factor is a monotonically decreasing function 
of the channel length. It also decreases drastically as Vsb is increased. Therefore, it 
has been modeled in this work as an exponentially decreasing function of Vsb- The 
expression for this model for the fudge factor 77 / is given in Eqn.(2.23). The value of 
the parameter rjfo appearing in the expression for rjf is a strong function of the effective 
channel length Leff and takes different values for different channel lengths, which are 
to be extracted from the subthreshold slope data. The coefficient k is almost constant 
with an approximate value of 0.75 for all channel lengths. The extracted values of 
the parameters 77/0 and k' for different channel length devices are given in Table 3.2. 
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Figure 3.5: Plot of Xg/Leff versus Leff for Vsb — 0-0 V for (a) Vds — 3.0 V and (b) 
Vds = 0.1 V. 
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The fudge factor rjf has been plotted as a function of Vsb in Figs.3.6(a) and (b) for 
the MOSFETs having effective channel lengths of 0.33 fim and 0.43 /rm respectively. 
In these figures, the fudge factor ’ Vialues calculated from its model [Eqn.(2.23)] are 
shown by continuous lines, whereas the values obtained from the measured values of 
the subthreshold slope from the experimental data are shown by the diamonds. An 
excellent match between the two is seen from the figures. 

Table 3.2 

The extracted values of the parzimeters used in the model of the fudge 
factor Tjf [Eqn.(2.23)] for the MOSFETs having effective channel lengths of 
0.33 /xm, 0.43 /xm, 0.63 ixta. zind 1.43 ^m. 


Channel Length (/xm) ; 

! VfO . 

k' (v)- 

1.43 

1.09 

0.75 

0.63 

^ 1.15 

0.735 

0.43 

1.54 

0.73 

0.33 

i 2.92 

0.75 


This model for the fudge factor is used for calculating the value of the minimum 
surface potential point Xo- The simulated results for the improved subthreshold slope 
model developed in this work[Eqn.(2.26)] are plotted as a function of Vsb in Figs.3.7(a) 
and (b) for Le// = 0.63 /xm and 0.33 /xm respectively for Vus = 0.1 V. Figures 3.8(a) 
and (b) show the same plots for Vos = 3 V. The experimentally measured data are 
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also shown in the figures for the sake of comparison. From these figures, it can clearly 
be seen that the improved model gives a much better fit to the measured data than 
the original model, which can be found by comparing these figures with Figs.3.3(c) and 
(d). Thus, our earlier assertion that the Xo point and thus the ideality factor rj are 
functions of applied Vds and Vsb is substantiated. 


3.2 The Subthreshold Current Characteristics 

In order to test the performance of the improved subthreshold slope model developed 
in this work, the Id-^gs characteristics are simulated for short channel MOSFETs 
in the subthreshold region. As described in the previous chapter, the subthreshold 
drain current model developed in this work[Eqn.(2.27)] uses the modified and improved 
subthreshold slope model based on the threshold voltage data given by the BSIM 
threshold voltage model[Eqn.(2.3)]. 

The performance of the subthreshold drain current model developed in this 
work is tested with the experimental data on short channel MOSFETs reported in 
literature [16,17]. The Id-Vcs characteristics of a MOSFET having an effective channel 
length of 0.075 /mi and an oxide thickness of 70 A have been given by Rittenhouse et 
al.[17] for Vsb = 0 V and for V^s values ranging from 0.05 V to 2.5 V in steps of 0.245 
V. 

In order to determine the values for the subthreshold slope 77, one needs to 
extract the values of the parameters 7, F, s, and S from the data of threshold voltage 
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Figure 3.6: The simulated plots for the fudge factor rjf versus Vsb for (a) Leff = 0.33 


ium and (b) Lgff = 0.43 fj,m. 


CENTRAL L4BR^ 

I 1 - T 


5k*. A Vi3^% 



Ideality Factor r| Ideality Factor r| 



0 12 3 4 

Source to Body Voltage (Volts) 

(a) 


I t I r - 1 I 1 



Source to Body Voltage Vsb (Volts) 

(b) 


Figure 3.7: The simulated plots of the ideality factor 77 obtained from the improved 
subthreshold slope model [Eqn.(2.26)] for = 0.1 V for (a) Te// = 0.63 72m and (b) 
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Figure 3.8; The simulated plots of the ideality factor rj obtained from the improved 
subthreshold slope model [Eqn.(2.26)] for V^s = 3.0 V for (a) Tg// = 0.63 /im and (b) 
L^ff = 0.33 fj.m. 
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variation with respect to the substrate potential Vsb- This requires the Id-Vgs char- 
acteristics for different values of V 55 . Unfortunately, these data are not presented by 
Rittenhouse et al.[17]. However, from Eqn.(2.26), it can be seen that the value of the 
parameter 77 can be extracted from the value of the parameter a, and in reality, the 
exact values of the parameters s and S are not needed. The values for the parame- 
ter a can be extracted from the threshold voltage versus V^s curve, which are found 
to be 0.435 for these cases. Similarly, the values of the parameters 7 and T can be 
extracted by simulating Eqn.(2.26) for one value of the measured subthreshold slope. 
For finding the value of x^/L, Eqn.(2.24) is used in which the value of the parame- 
ter k' has been extracted earlier. The subthreshold slope model results are plotted in 
Fig.3.9, along with the experimentally measured values for the sake of comparison. An 
excellent match is seen from the figure, which stakes the claim that the subthreshold 
slope model developed in this work is very accurate even for the sub- 0.1 micron channel 
length MOSFETs. This model of the subthreshold slope has been used in the drain 
current model. 

The simulated Id-Vgs characteristics in the subthreshold region using the im- 
proved subthreshold slope model for this device are shown in Fig.3.10, where the exper- 
imentally obtained results[17] are also shown for the sake of comparison. From Fig.3.10, 
it is seen that the simulated results obtained from the subthreshold drain current model 
developed in this work show a good match with the experimental data, particularly 
in the deep subthreshold region. In the moderate inversion region, the performance of 
our model deteriorates. This is a problem of the global parameter extraction approach 
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Figure 3.9: Comparison of the subthreshold slope calculated from our improved sub- 
threshold slope model [Eqn.(2.26)] with the subthreshold slope obtained from the ex- 
perimentally measured data reported by Rittenhouse et al.[17] for a 0.075 ^m channel 
length device. 

as has been elaborated by Deshpande and Dutta[10]. 

The simulated results of our subthreshold slope model and the Id-^gs charac- 
teristics have been compared with the experimental data reported on another sub-0.1 
micron MOSFET, where the effective channel length was 0.09 /rm, and the oxide thick- 
ness was 35 A [16]. The Id-Vqs characteristics have been given for Vbs 0.05 V and 
1.55 V by Mii et al.[16]. The values of the model parameters 7, F, s, and S are extracted 
by the method similar to the one applied for the 0.075 /xm channel length device. The 
subthreshold slope model results are plotted in Fig.3.11 along with the experimentally 
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Figure 3.10: Comparison of the Id-Vos characteristics calculated from our model with 
the experimental data reported by Rittenhouse et al.[17] for a 0.075 jjm channel length 
device. 

obtained values[16] for the sake of comparison. It is observed that the match between 
the two is quite good. 

The simulated Id-Vgs characteristics in the subthreshold region for the 0.09 
/iin channel length MOSFET reported by Mii et al.[16] are shown in Fig.3.12, where 
the experimentally obtained results are also shown for comparison. This figure shows 
a good match between the simulated results and the experimentally obtained results. 
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Figure 3.11: Comparison of the subthreshold slope calculated from our improved sub- 
threshold slope model [Eqn.(2.26)] with the subthreshold slope obtained from the ex- 
perimental data reported by Mii et al.[16] for a 0.09 ^m channel length device. 


3.3 Results of the Complete Drain Current Model 

The complete model for the drain current, valid for all regions of operation, has been 
described in the previous chapter. The threshold voltage model and the subthreshold 
slope model developed in this work are used in the complete model given by Deshpande 
and Dutta[10]. The complete model, described in detail in Section 2.8, has been simu- 
lated for the 0.075 ^m and 0.09 ixm channel length MOSFETs reported by Rittenhouse 
et al.[17] and Mii et al.[16] respectively. The complete Id-Vgs characteristics simulated 
from the model are plotted for the 0.075 channel length MOSFET in Fig.3.13 along 
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with the experimentally obtained values[l7]. Similarly, same characteristics are plotted 
in Fig. 3. 14 for the 0.09 ^m channel length MOSFET reported by Mii et al.[16]. The 
results show a very good match. 

The simulated Jp-Vbs characteristics are plotted in Fig.3.15 for different values 
of the applied gate voltage Vos for the 0.075 channel length MOSFET, along with 
the experimentally obtained results reported by Rittenhouse et al.[17]. 



Figure 3.12; Comparison of the Id-Vgs characteristics in the subthreshold region cal- 
culated from our model with the experimental data reported by Mii et al.[16] for a 0.09 
jum channel length device. 
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Figure 3.13: Comparison of the Id-Vgs characteristics calculated from our model with 
the experimental data reported by Rittenhouse et al.llT) for a 0.075 ,m. channel length 

device. 
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Figure 3.14; Comparison of the Id-Vqs characteristics calculated from our model with 
the experimental data reported by Mii et al.[16] for a 0.09 [xui channel length device. 
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Figure 3.15: Comparison of the Id-Vds characteristics calculated from our model with 
the experimental data reported by Rittenhouse et al.[17] for a 0.075 nm channel length 
device. 
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Chapter 4 


Summary and Conclusion 


The emergence of the analog integrated circuits and mixed-signal VLSI chips demands 
very accurate and elaborate MOSFET models suitable for circuit simulation. The 
MOSFET models which satisfy the requirements of digital circuit simulation are cur- 
rently no longer useful for analog circuit simulation. Therefore, MOSFET mode lin g 
suitable for analog circuit simulation has become one of the most challenging fields in 
the VLSI area. Also, very high packing density achieved by the reduction of the size of 
an individual transistor in the IC chip has become possible by the improvements in the 
photolithography techniques. This high packing density created high power dissipa- 
tion problem. Thus, low operating current and voltage designs are becoming popular 
because of their inherent advantages in reducing the power dissipation, which needs to 
be kept within limit. The low voltage operation pushes the MOSFET into its weak 
inversion region. In addition to this, quite a few battery operated applications use the 
MOSFETs in their subthreshold region. Hence, it has become imperative to accurately 
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model the MOSFET in its subthreshold region. 


4.1 Summary of the Work 

The work done is summarized in this section. Initially, the subthreshold slope model 
proposed by Kang et al.[9] has been simulated in this work. For this simulation, the 
various parameters used in the model are extracted from the threshold voltage data. 
The simulation of this subthreshold slope model is done for MOSFETs having effective 
channel lengths of 0.33 iim, 0.43 urn, 0.63 //m, and 1.43 ;um. These simulated results 
have been compared with the experimentally obtained results of the subthreshold slope 
from the data reported by Lee and Min[13] for these device. It has been observed from 
this simulation that the simulated results of the subthreshold slope overestimate its 
measured values, more so for the MOSFETs having relatively shorter channel lengths 
(e.g., for Leff = 0.43 /^m, 0.63 /xm, and 0.33 /xm). This is due to the improper modeling 
of the drain voltage dependence of the subthreshold slope by Kang et al.[9]. 

This subthreshold slope model is based on the capacitor equivalent circuit of 
the MOSFET in the subthreshold region. The subthreshold slope equation of this 
model is analytically expressed in terms of the parameters of the BSIM threshold 
voltage model[6] by carefully interpreting it with the capacitor equivalent circuit of 
the MOSFET in the subthreshold region. In this model, it has been assumed that 
the distance of a point on the channel from the source where the surface potential is 
minimum converges rapidly to a value equal to L^ff/Z as the drain voltage is increased. 
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Thus, the dram voltage dependence of this minimum surface potential point has been 
neglected by Kang et al.[9]. Hence, the subthreshold slope expression obtained in their 
work comes out to be independent of the drain voltage. However, experimental data on 
the short channel devices reported in the literature show that the subthreshold slope is 
a strong function of the drain voltage. Thus, this model has a major flaw, which makes 
it highly unsuitable for the simulation of short channel MOSFETs in their subthreshold 
region of operation. In this work, this model has been modified in order to develop 
a more accurate subthreshold slope model which takes into account the effects of the 
drain voltage, as well as that of the effective channel length, and the body voltage on 
the subthreshold slope. 

In order to develop this improved model, the distance of the minimum surface 
potential point Xo is evaluated for various values of Le// and Vbs from the analytical 
expression given by Liu et al.[14], after making some modifications, as detailed earlier 
in Chapter 2. This model of Xo has been used in the original subthreshold slope model 
developed by Kang et al.[9]. Thus, a new and improved subthreshold slope model is 
obtained in this work. The values of the parameters 7, T, s, and S used in the improved 
subthreshold slope model have to be extracted from the threshold voltage data. These 
parameters are used to model the various short channel effects. In addition to these 
four parameters, this model has one more fitting parameter rjj, whose value is to be 
extracted from the measurement data for short channel MOSFETs. The simulation of 
this improved subthreshold slope model is done for MOSFETs having effective channel 
lengths ranging from 0.075 pm to 1.43 pm, and the simulated results are observed to 
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be matching quite closely with the measured values of the subthreshold slope. This 
improved subthreshold slope model is based on a true physical foundation of device 
physics. It has been made accountable for the DIBL effect along with the dependence of 
the subthreshold slope on the various parameters such as the effective channel length, 
the body voltage, and the drain voltage. 

The drain current model of Deshpande and Dutta[10] is modified in this work 
by replacing the subthreshold slope model and the threshold voltage model used in 
it by the improved subthreshold slope model developed in this work and the BSIM 
threshold voltage model respectively, in order to arrive at a new drain current model. 
This new drain current model is simulated for MOSFETs having effective channel 
lengths of 0.075 /im and 0.09 /um reported by Rittenhouse et al.[17] and Mii et al.[16] 
respectively. The results obtained from this simulation are observed to be matching 
very well with the experimentally obtained results. 


4.2 Conclusion 

1. The subthreshold slope model developed by Kang et al.[9] is derived after as- 
suming that the value of Xo (i-6., the distance of a point on the channel from the 
source where the surface potential is minimum) is equal to Le///3, independent 
of the drain voltage and the effective channel length. However, this is not true. 
The value of Xo is observed to be strongly dependent on the various parameters 
such as the drain voltage, the body voltage, and the effective channel length. 
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Therefore, in order to model the effect of the drain voltage on the subthreshold 
slope, it is important to model Xo in terms of these parameters. 


2. In this work, in the calculation for the value of Xo, we have used a fitting parame- 
tei, referred to as the fudge factor rjf. This has been done in order to compensate 
for the invalid assumption of a constant depletion layer thickness along the chan- 
nel in the subthreshold region. It has been observed that the value of the fudge 
factor decreases as the effective channel length is increased, and it finally con- 
verges to a value equal to unity for MOSFETs having longer channel lengths. 
Also, the value of the fudge factor decreases exponentially with an increase in 
the body voltage Vsb- Though the fudge factor is modeled as a fitting parameter 
in this work, however, its behavior with respect to the effective channel length 
and the body voltage implies that it represents a true physical phenomenon due 
to the non-uniform depletion layer thickness along the channel of the MOSFET. 
Thus, this parameter rjf has a well justified physical meaning. 


4,3 Future Scope of Improvement 

As explained in the previous section, the fudge factor rjf has a physical meaning. 
However, in this work its value has to be extracted from the measurement data, and the 
exact dependence of this fudge factor value on the substrate doping, effective channel 
length and oxide thickness can not be defined at this moment. The sub threshold model 
developed in this work can be improved by analytically expressing this fudge factor in 
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terms of the terminal voltages and the device specifications such 
length, the oxide thickness, and the substrate doping. 


as the effective channel 
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